Redundancy of classical and quantum correlations during decoherence 
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We analyze the time dependence of entanglement and total correlations between a system and 
fractions of its environment in the course of decoherence. For the quantum Brownian motion model 
we show that the entanglement and total correlations have rather different dependence on the size 
of the environmental fraction. Redundancy manifests differently in both types of correlations and 
can be related with induced-classicality. To study this we introduce a new measure of redundancy 
and compare it with the existing one. 
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I. INTRODUCTION 

Decoherence is a physical process taking place when a 
system interacts with its environment. This process is 
crucial to understand the origin of the classical domain 
from a fundamentally quantum substrate [T]. Absence of 
macroscopic superpositions is one of the essential features 
of the classical realm which is explained by decoherence 
as a consequence of an effective super-selection rule that 
prevents the stable existence of the vast majority of states 
in the Hilbert space. In recent years, important devel- 
opments enabled us to understand the origin of another 
defining property of classical systems: the fact that they 
exist in an objective state, a state of their own. Zurek 
and co-workers made this idea precise by noticing that 
the emergence of classicality is connected with the ex- 
istence of redundant records of the state of the system 
imprinted in the environment. Redundancy turns out to 
be the key notion to define objectivity and classicality: A 
system behaves classically when different fractions of any 
environment correlate with it in the same way (i.e. when 
correlations become redundant). Such consensus about 
the state of the system characterizes the classical realm 
and enables a definition of objectivity. In their pioneer- 
ing work, Zurek and co-workers [H |21 SI E E] analyzed 
the emergence of redundancy in the total correlations be- 
tween the system and the environment which can be mea- 
sured by using information theoretic tools such as mutual 
information. Here, we will examine how does redundancy 
manifests in purely quantum correlations (entanglement) 
established between the system and fractions of the en- 
vironment in the course of decoherence. We will present 
a definition of redundancy based on the entanglement 
and show that the way in which redundancy manifests in 
total and quantum correlations is rather different. 

We will focus on the paradigmatic quantum Brownian 
motion (QBM) model H] where a particle S interacts 
with an environment £ that can be split into a fraction 
£f and its complement £i-f (/ parameterizes the size 
of the fraction £f with respect to £ ranging between 
and 1). Total correlations between S and £f can be mea- 
sured with the mutual information T{S,£f) (computed 
from the reduced state obtained after tracing out the 
complement fi_/). On the other hand, quantum corre- 
lations can be quantified with an entanglement measure 



between S and £f, denoted as E'^^\ Here we will study 
both I{S,£f) and E'^^^ as a function of time and /, for 
random choices of the components of £f out of the total 
£. 

The paper is organized as follows: In Sec. |ll]we present 
the model and describe the temporal evolution of the cor- 
relations between the system and the different portions of 
the environment. We consider two types of splittings of 
the environment in different fractions. First we analyze 
a splitting into bands oscillators (characterized by their 
frequencies) and then we analyze grouping such bands 
into random fractions (characterized by their size). In 
Sec. |III| we derive analytical results for the mutual infor- 
mation and entanglement in the non-dissipative regime. 
In Sec. IV we study the evolution of redundancy of to- 
tal correlations and quantum correlations. In Sec. [V| we 
summarize our results. 



II. DYNAMICS OF CORRELATIONS IN 
QUANTUM BROWNIAN MOTION 

In the QBM model a central particle S is coupled 
to an environment £ formed by harmonic oscillators q„ 
{n = l,..,iV). The total Hamiltonian is H — Hg ' 



Hs£ + Hs where Hg = /2m + TOCj|a;2/2, Hg = 
E„ (2^ + and Hs£ = a;E„c„g„ (we use 

h ^ 1 throughout). The spectral density J{uj) ~ 
c^(5(ti; — LOn) /2nnnijJn determines the effect of £ on 
S . We consider a realistic family of environments where 
J(w) = 277170^ 9{K — lo)/ti that includes Ohmic 

(n = 1), super-Ohmic, (n > 1) and sub-Ohmic (n < 1) 
members (A is a high frequency cutoff, 70 a coupling 
strength, and 6{x) the Heaviside step function). We 
assume that the initial state of £ is the ground state 
of and that system-environment correlations vanish 
initially, when the interaction is switch on correlations 
are created between system and environment [3] . Quan- 
tum mutual information (MI) X{S,£f) (TUl |TT] is de- 
fined from the von Neumann entropies of the correspond- 
ing states: 1(5, £f) = n{S) + n{£f) ~ 'H{S,£f) (here 
Ti-iS) — —Tr[pslnps], etc.). The reduced state P{s,£f} 
is obtained after tracing out the complementary environ- 
ment £i-f. We restrict to initial Gaussian states which, 
due to linearity, remain Gaussian for all times and can 
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be described efficiently by a covariance matrix. In such 
case analytic expressions for T{S, £f) can be obtained. In 
order to quantify the quantum correlations, we will use 
the multipartite logarithmic negativity [12^ as a measure 
of entanglement between S and £f. This is defined as the 
maximum between zero and — X^i iy <i/2 ^'-'SC^^'i); where 
i>i are the symplectic eigenvalues of the partially trans- 
posed covariance matrix corresponding to the Gaussian 
state P{s,£f}- 

We present first numerical results and later we show 
analytic expressions that reproduce them with high ac- 
curacy. To do numerics we use a discrete version of the 
model [H [T3] with 600 modes (we checked that our re- 
sults are stable as this number is increased) and coupling 
strengths proportional to the spectral density. We set 
all masses to unity and take 70 = 0.1, the renormalized 
frequency ils = 3, the cutoff frequency A = 20 for dissi- 
pative environments (n = 1, 1/2) and we use a high cut- 
off, A = 300, for the non-dissipative one (n = 3). The ini- 
tial state of the system is a pure Gaussian squeezed state 
with a squeezing parameter r — In^mflgAx/ Ap) = —5 
(unless it is specified) and AxAp — 1/2. 

Correlations with bands of oscillators: Mutual infor- 
mation and entanglement between the system and envi- 
ronmental bands of frequency w are shown in Fig. [T] MI 
develops fast and for short times concentrates at high fre- 
quencies. For dissipative environments (sub-Ohmic case, 
n = 1/2) at later times the resonant band is dominant. 
After a few relaxation times, the resonant peak is washed 
out and all oscillators in the environment tend to play 
a similar role. The entanglement has a similar behav- 
ior except for the super- Ohmic case. The resonant band 
of the environment is dominant for entanglement and is 
also washed out by dissipation (contrary to the case of 
mutual information, entanglement with non-resonant os- 
cillators seems to be always negligible). One of our find- 
ings is that the super-Ohmic environment displays re- 
coherence effects which manifest as oscillatory behavior 
both for the mutual information and entanglement. In 
this case, entanglement only develops with resonant os- 
cillators. Thus, the super-Ohmic environment does not 
induce truly dissipative effects but is responsible for a 
renormalization (dressing) of the system as well as for 
coherent oscillations which will also be visible in our 
forthcoming studies. It is worth pointing out that the 
existence of a perfect reversion for the super-Ohmic en- 
vironment only takes place for some initial sates of the 
system. If instead of the above state we choose one with a 
large momentum squeezing complete recoherence is lost 
(see below). This was already noticed in [14| and can 
be attributed to a process taking place at ultra-short 
timescales, as will be evident below. 

Correlations with fractions of the environment: We use 
"partial information plots" (Pl-Plots) and "partial en- 
tanglement plots" (PE-Plots) to study the evolution of 
correlations between S and fractions £f of variable size /. 
They are plots of the average mutual information (aver- 
age of MI over random choices of fractions 5/ of the same 
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FIG. 1: Mutual Information (left) and logarithmic negativity 
(right) between a Brownian particle and bands of oscillators 
of the environment of frequency ui. The initial state of the 
system is squeezed in position (r = —5). Plots at the top cor- 
respond to a sub-Ohmic environment while the ones at the 
bottom correspond to the super-Ohmic case. Ohmic environ- 
ment behaves like sub-Ohmic one. 



size) and average entanglement between the S and £f, 
respectively. The evolution of the Pl-Plots for different 
spectral densities are shown in Fig. [2] (we verified that for 
the Ohmic and sub-Ohmic cases the results are similar), 
where for simplicity we subtract the system's entropy and 
plot T{S,£f) — 7i(5). As the total state is pure, mutual 
information is l{S,£f) — H{S) + H{S, £i-f) — H{S,£f), 
i.e. it is symmetric around / = 1/2 [3]. This means that 
by adding the mutual information between the system 
and two complementary fractions of the environment, 
one always obtains the maximal available information 
I{S,£) = 2n{S) (i.e. I{S,£f) +I{S,£i-f) = 2n{S)). 
Redundancy clearly manifests in the Pl-Plots: There is 
a sharp growth of the mutual information for small / 
followed by a plateau that indicates that mutual infor- 
mation becomes almost independent of the fraction /. 




FIG. 2: Partial Information Plots (Pl-Plots) for a dissipative 
(sub-Ohmic) environment (left) and non-dissipative (super- 
Ohmic) one (right). Initial state is squeezed in momentum 
with r = —5. The super-Ohmic case displays clear recoher- 
ence events for which the system and the environment loose 
their correlations. 
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As Pl-Plots are not completely flat around / = 1/2 
there is also some non-redundant information Xj^ir [4 
that can be defined as the slope of at / = 1/2 

(see below). Our results show that for dissipative envi- 
ronments (both Ohmic and sub-Ohmic) Iatt?, is dramati- 
cally reduced, much faster than the total available infor- 
mation. The behavior of the super-Ohmic environment 
is rather different: redundancy develops but this effect 
is reversible. Decoherence is followed by recoherence and 
de-correlation. Thus, the super-Ohmic environment does 
not induce truly dissipative effects but is responsible for a 
renormalization (dressing) of the system as well as for co- 
herent oscillations which will also be visible in our forth- 
coming studies. Perfect reversion for the super-Ohmic 
environment only takes place for some initial states of 
the system. If instead of the above state we chose one 
with a large momentum squeezing complete recoherence 
is lost. This was already noticed in [14] and is due to 
transitory effects due to the instantaneous switching of 



the interaction. 

Partial entanglement plots (PE-Plots) in Fig. [sjshow 
that for small (and moderate) values of /, entanglement 
scales almost linearly with / and stays very small. For 
larger fractions the rapid growth of entanglement in- 
dicates that it is established between the system and 
the environment " as a whole" . The effect of dissipa- 
tion is displayed in the PE-Plots: Dissipation reduces 
the maximum attainable entanglement and flattens the 
initial slope of the PE-Plot, suppressing even more the 
quantum correlations with small fractions of the envi- 
ronment. Super-Ohmic PE-Plot is dramatically different 
from the one corresponding to the dissipative case. In 
this non-dissipative environment revivals are evident in 
the PE-Plot. Redundancy of quantum correlations is re- 
lated with the flatness of the PE-Plot as it indicates that 
many small environmental pieces become entangled as a 
whole with the system S . A way to quantify this notion 
of redundancy will be presented below. 
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FIG. 3: Partial Entanglement Plot (PE-Plot) for sub-Ohmic (left) and super-Ohmic environments (center and right). The 
initial state of the system is squeezed in position with r = —5 except for the right plot where the squeezing is in momentum 
r = 5. In that case recoherence is suppressed. Entanglement grows very slowly as a function of the fraction size /, which is 
evidence of the branch structure of the total state. 
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III. ANALYTICAL RESULTS 

A simple model enables us to reproduce and under- 
stand the numerical results: We assume that S is mas- 
sive and under-damped, justifying a Born-Oppenheimer 
approximation [5]. Under this approximation, the evo- 
lution of the joint state oi S ~ £ has a branch structure 
that can be written as: 



I*; 



S£ 



= / il^{x,t)\x)s®\Mi))£i ■ ■ M'>Px{t))s„dx, (1) 



where for each environmental oscillator, the state 
\i^x{t))£„ depends parametrically on the system's posi- 
tion X. Its time evolution can be found by noticing that 
each oscillator feels an external force F„(t) — Cnx{t). 
Here, x(t) parametrically depends on x and has a rather 
simple form for large squeezings: For positive (nega- 



tive) r a trajectory starts from a;(0) = x {x{0) — 0) 
with velocity i(0) = (i(0) = fls^)- Thus, in gen- 
eral, x{t) = 9{r)x cos{^lst) + 6{—r)xsm{^lst). Us- 
ing this, the fate of initial squeezed Gaussian states 
can be found by computing all elements of the co- 
variance matrix of the full state. The entanglement 
Eif) can be analytically computed and turns out to 
depend upon the function S-^f\t) — J2ne£f'^n — 
T,n€£f 4 {Ko-Ut) + al{t)) /4to„w„, where 



{C0s{nst) — COs{Wnt)) . 
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On average, this satisfies d^^f\t) = f d{t) and 

E<^f'\t) = -^lii^l+4d{t)Sx^{l + 3f) (2) 

+ dit)5x^{l + 3f) ] - ( ^^V}- 



,1 + 3/ 



1 + 3/ 



where Sx is defined in terms of the absolute value of 
the squeezing parameter |r| — hi(Sx / Sxq) being Sxq the 
ground state variance Sxq = {l/2mflsY^^ . This re- 
sult reproduces with high accuracy the numerical data 
before dissipation sets in, see Fig. |4] and replacing 
/ — ^ dn{t)/d(t) we obtain the entanglement with the 
band ?«„. Before analyzing this, we present a similar 
result for the mutual information [5], that turns out to 
be 

I{S,£f) = h{x{l)) + hixif)) - h{x{l - /)). (3) 

Here h{x) is the entropy function h(x) = (x+1/2) ln(x+ 
1/2) — (x— 1/2) ln(x— 1/2) and its argument xif) satisfies 



X(/) 



2/ d{t)Sx'^ 



(4) 



The three terms in the equation ([3| are, respectively, the 
entropies for the system, the environment and the joint 
S — £f ensemble (notice that xif = 1) is the symplectic 
area of S). It is worth noticing that this simple model 
enables us to explain the peculiar super-Ohmic case (with 
high cutoff) for which d{t) is 



d{t) 



2tt 



{9{-r) sin^ m + e{r){l + cos^ nt)) . (5) 



We can use this in the above formulae for E'^^'^ and see 
that a dramatically different behavior for positive and 
negative squeezing arises, fully reproducing the numer- 
ical data. For initial states delocalized in position the 
sudden switch of the interaction produces a kick for each 
environmental oscillator whose effect cannot be undone 
later (the environmental oscillators remain out of phase 
and no recoherence is possible). On the other hand, for 
initial states delocalized in momentum the initial kick is 
negligible and the transient effects are minimized. An 
equivalent behavior is observed if the interaction is adi- 
abatically switched on. The existence of recoherence is 
not a novel result |14j , but is remarkable that it shows up 
in such a clear way for the entanglement and the mutual 
information. 

£^(^^ and 1(5, £f) have similarities but also remarkable 
differences. E^^^ grows slowly for small / and fast for 
/ w 1. For d{t)Sx^ 3> 1, we have 



« - In 
2 



(1 + 3/)= 



(l-/)(l + 3/)2 + 2//d(<)(5a;2 



(6) 



Numerical results are compared with the analytic esti- 
mate of E^^'^ in Fig. |4| Notably, in the limit of large 



squeezing and / < 1 the entanglement E^^'' becomes in- 
dependent of the environmental spectral density that en- 
ters the above expressions through the function d{t). In 
fact, in such case we have 



1, {1 + 31 



(7) 



a similar behavior was observed for Gaussian 
Greenberger-Horne-Zeilinger (GHZ)-type states |15j . 
In the presence of dissipation i?'^^ is reduced so the 
above universal expression is really an upper bound 
for E^^\ A remarkably similar result arises for the 
mutual information. The analysis in this case is slightly 
different since MI depends on the squeezing through 
the combinations fd{t)Sx^ and (1 — f)d{t)Sx^. Thus, 
both for small and large fractions the first derivative 
of X(5,£/) is singular (due to the singularity in h{x))- 
In the limit d{t)5x^ > 1, n{S) « hi{2e^ d{t) 5 x"^) / 2 and 
mutual information behaves as 



i{s,£o<}<i)^n{s) 



1 



In 



/ 



1-/ 



(8) 



X(5, £f) grows fast for small / and approaches a plateau 
with a slope related to the so-called non-redundant in- 
formation Inb.- This also turns out to be universal since 
Ink = dI{S,£f)/df\f^,/2 « 2, see Fig. |4) 
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FIG. 4: Top: Normalized PE-Plot and Pl-Plot for a fixed 
time. Analytical estimates (solid lines) given by expressions 
in equations ([2| and (|3| show a remarkable agreement with 
numerical data (dots). The relevant fractions /e and // to 
compute information and entanglement redundancies Re and 
Ri are shown in the figure. Bottom: Re and Ri as a function 
of time are shown for a dissipative environment (n — 1/2), 
where the relaxation timescale is treiax ~ 3 (with a deficit 
Se = 25i = 0.2). 
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IV. REDUNDANCY OF CORRELATIONS 

Our results for E^^") enable us to define a new measure 
quantifying the redundancy of quantum correlations. We 
define the entanglement redundancy, Re = l/Zs, as the 
number of environmental fractions of size fs such that 
by ignoring one of them we induce a decay of the aver- 
age entanglement to a fraction 6e of the maximal one. 
Roughly speaking, a large entanglement redundancy is 
obtained when by ignoring a small fraction of the envi- 
ronment one looses a large fraction of the total entan- 
glement. This is the case if /e turns out to be small 
when 5e small. GHZ states have large entanglement 
redundancy, since by tracing over one of its constituents 
the entanglement is lost. Indeed, large entanglement re- 
dundancy characterizes the fragility of the multipartite 
entanglement. Computing Re from the PE-Plots is sim- 
ple: As shown in Fig. jl] we should find the smallest 
fraction Je such that E^^'^^ = SeE^^='^\ Since Re is 
the number of such fractions in the environment £, we 
have Re = I/Zb- As E^^^ is a monotonous function of 
f, a large redundancy is expected for branch states eq. 
M, in that case E^^'' would grow very slowly with / for 
small and moderate values of /. Remarkably, our ana- 
lytic results before dissipation enable us to obtain a good 
estimate for Re'- 

Re^^^c''-"^''^"' ^A{tr-, (9) 

JE 

where A{t) = d{t)5x^ / SqX^ is the symplectic area of 
the system's state (in the limit of large squeezings) and 
^(/=i) ~ ln(32(i(i)(5a:^)/2 is the maximal entanglement 
(which is proportional to the entropy of the system). 

It is interesting to notice that entanglement redun- 
dancy Re can be related with information redundancy 
Ri, defined in [5 as the number of environmental pieces, 
Ri = 1/fi, that carry a fraction (1 — (5/) of the available 
classical information H{S) (see Fig. [4|. In Fig. [3] we 
can see that both redundancy measures. Re and Rj, be- 
have in a similar way. Notably, both Re and Ri grow for 
times much longer than the relaxation timescale. This is 
a consequence of the fact that when dissipation becomes 
effective the quantum correlations with fractions of the 
environment are almost erased. Consequently, PE-Plots 
and Pl-Plots flatten well before equilibration. Finally, in 
the asymptotic limit, redundancy almost disappears as 



do all types of correlations when the system approaches 
equilibrium. Our results, also show that the amount of 
non-redundant information is intimately related with the 
entanglement created between S and fractions of the en- 
vironment. Besides, we can show that both redundancies 
numerically coincide when the deficits Se and Sj are such 
that 6e ~ 5in{S)/E^f=^^ i;(/=i/2)/i;(/=i). Thus, in 
the limit of large squeezing both deficits become identical 
(this is the case because i?(/=i/2) is bounded from above 
asi;(/=i/2) <lny5). 

V. CONCLUSIONS 

In this paper we studied the evolution of quantum cor- 
relations (characterized by the entanglement) and total 
correlations (measured by the MI) between a system and 
its environment. First we analyzed the evolution of such 
correlations as a function of the band's frequency. For 
this case we showed that for dissipative environments 
the resonant bands are always the ones that correlate the 
most with the system. On the other hand, for the non- 
dissipative environments the high frequency bands dom- 
inate the evolution inducing short time effects that can 
leave a long time footprint for some initial states. Other 
initial states lead to revivals of the correlations. Our ana- 
lytic and numerical study of the evolution of correlations 
for QBM shows that MI and entanglement do become re- 
dundant at a time-scale which is connected with the one 
of decoherence (indeed, the function d{t) appearing above 
is the same that characterizes the exponential suppres- 
sion of the quantum interference effects in QBM). This is 
seen in the Pl-Plots and PE-Plots. Our results also show 
the special nature of the state that dynamically evolve for 
the QBM model that has a branch structure compatible 
with large entanglement redundancy and truly multipar- 
tite entanglement between S and £. Dissipation rapidly 
erases quantum correlations with fractions of the envi- 
ronment and flattens both Pl-Plots and PE-Plots. Thus, 
redundancy is dramatically increased by the action of 
dissipation. 
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